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ABSTRACT 

Coronal mass ejections (CMEs) are large-scale ejections of plasma and mag- 
netic field from the solar corona, which propagate through interplanetary space 
at velocities of ~100-2500 km s -1 . Although plane-of-sky coronagraph measure- 
ments have provided some insight into their kinematics near the Sun (<32 R ), 
it is still unclear what forces govern their evolution during both their early ac- 
celeration and later propagation. Here, we use the dual perspectives of the Solar 
TErrestrial RElations Observatory (STEREO) spacecrafts to derive the three- 
dimensional kinematics of CMEs over a range of heliocentric distances (~2- 
250 R ). We find evidence for solar wind (SW) drag-forces acting in interplan- 
etary space, with a fast CME decelerated and a slow CME accelerated towards 
typical SW velocities. We also find that the fast CME showed linear (5 = 1) 
dependence on the velocity difference between the CME and the SW, while the 
slow CME showed a quadratic (5 = 2) dependence. The differing forms of drag 
for the two CMEs indicate the forces and thus mechanism responsible for there 
acceleration may be different. 

Subject headings: Sun: coronal mass ejections (CMEs) — Sun: activity — solar- 
terrestrial relations 



1. Introduction 

Massive eruptions of plasma and magnetic field which travel from the Sun through the 
Heliosphere are known as coronal mass ejections (CMEs). CMEs can have masses up to 



10 16 g (Vourlidas et al. 2002), propagate at velocities of up to 2500 km s 1 (Gopalswamy 
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2004) close to the Sun, while at 1 AU velocities tend to be closer to that of the solar 
wind (SW; |Gopalswamy 2007). Although CMEs have been the subject of study for nearly 
40 years, a number of fundamental questions regarding their acceleration and propagation 
remain unanswered. One such question, what forces govern the propagation of CMEs in the 
Heliosphere has been especially difficult to tackle. This is mainly due to a lack of the three 
dimensional (3D) observations of CMEs in the inner Heliosphere. 

The kinematic evolution of CMEs can be broken into three phases; initiation, accel- 
eration, and propagation (Zhang et al. 2001). During the propagation phase, the initial 



acceleration has ceased and the CME motion is dominated by the interaction between the 
SW and the CME. The "snow plough" , aerodynamic drag, and flux-rope models all aim to 



explain the motion of CMEs in the SW (Tappin 2006 Borgazzi et al. 2009; Vrsnak et al. 



2010 


Cargill 


2004; 


Chen 


1996) 



dominated regime may be written: 

dv 
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where Co is the drag coefficient, p sw is the solar wind density, A cme is the CME area and v sw 
is the solar wind velocity, and M is the CME mass. We use a parametric drag model similar 



to that of Vrsnak & Gopalswamy (2002) with the added parameter 5, which determines if 
the drag is quadratic or linear. This parametric form collapses the complex dependences of 
the CME area (A cme ) and the solar wind density (p sw ) into a power-law which depends on 
heliospheric distance R. Eq. [I] can thus be written 

dv 
li 

where a, (3, and 5 are constants. 
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Before the launch of the Solar TErrestrial RElations Observatory (STEREO; Kaiser 



et al. 2008 ) mission, synoptic white-light CME observations were limited to 32 R using 



Large Angle Spectrometric Coronagraph (LASCO; Brueckner et al. 1995), while the Solar 



Mass Ejection Imager (SMEI; Jackson et al. 2004, Howard et al. 2006) sometimes tracked 
CMEs to Earth (~215R Q ). In radio observations, fast CMEs which drove shocks could be 



tracked to Earth (Reiner et al. 2007). Interplanetary Scintillation (IPS) observation pro- 
vided density and velocity measurements for both CMEs and the SW from 50 R to beyond 



1AU and using tomographic techniques can give 3D information (Manoharan 2006, 2010). 
CMEs are also observed in in- situ measurements with WIND and ACE at LI (~1 AU), and 
occasionally CMEs can be tracked up to very large distances of up to 5 AU using additional 



spacecraft (Tappin 2006). Numerical modelling has been used to study CME propagation 



with numerous approaches such as, ID Hydro simulations, 2.5D MHD simulations and full 
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3D MHD simulations ( 


Gonzalez-Esparza et al. 


2003 


Cargill et al. 


1996 


Cargill & Schmidt 


2002 


Odstrcil & Pizzo 


1999 


Odstrcil et al. 


2004 


Smith et al. 


2009 


Falkenberg et al. 


2010 


)• 



Statistical studies comparing in-situ with white light observations indicate a trend of 



CME velocity converging towards the SW velocity as they propagate to 1 AU (Gopalswamy 



2007). Other studies, based on white light observations have indicated that aerodynamic 



drag of some form may explain this trend (Vrsnak 2001 Shanmugaraju et al. 2009). Radio 



observations suggest that a linear form of aerodynamic drag is most appropriate for fast 
CMEs ( |Reiner et al. 2003). Tappin (2006) showed that acceleration can continue far out 
(5 AU) into the Heliosphere. However, these studies are subject to the difficulties associated 
with the observations they are based on. For example white light observations were limited 
to single, narrow, fixed, view-points meaning only observation of the inner Heliosphere could 



be made and even these were subject to projection effects (Howard et al. 2008b). Also, 
linking features in imaging and in-situ observations is complex and can be ambiguous, a 
problem exacerbated during periods of high activity. In the case of numerical simulations 
their complexity can make it hard to extract which effects are the most important, possibly 
obscuring the important underlying physics. 

The unique STEREO mission consists of two nearly-identical spacecraft in heliocentric 
orbits, STEREO-B(ehind) and STEREO-A(head) which separate from the Sun-Earth line at 
22.5° per year. Each spacecraft carries the Sun Earth Connection Coronal and Heliospheric 



Investigation (SECCHI: Howard et al. 2008a) suite, which images the inner Heliosphere 



from the Sun's surface to beyond 1 AU. Using STEREO observations, a number of papers 
have been published which extract 3D information and study CMEs at extended heliocentric 



distances, over-coming some of the difficulties outlined above. Davis et al. (2009) identified a 
CME in HI1 and HI2, using a constant velocity assumption ( Sheeley et al.|2008 ) they derived 
the speed and trajectory of the CME. The predicated arrival time, based on the speed derived, 



agreed with the in-situ observations. Wood et al. (2009) used the "Point-P" and "Fixed- 



methods to derive the height, speed, and direction from elongation measurements out to 
distances of ~12OR . A recent paper by Liu et al. (2010) tracked a CME to ~15OR in 3D 



using J-maps from both spacecraft to triangulate the CMEs position in 3D. On the other 



hand Maloney et al. (2009) tracked the trajectory of CME apexes in 3D using triangulation, 



some as far as 240 R . Byrne et al. (2010) developed a new reconstruction method with 
allowed the entire CME front to be reconstructed. They found evidence for CME deflection, 
expansion and acceleration low down (<7R ) followed by a solar wind drag interaction. 
For a review of some of the different 3D reconstruction methods which have been applied to 



STEREO CME observations see Mierla et al. (2010). 



In this paper, we use triangulation to localise CME apexes in 3D. From this, we derive 
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the CME apex trajectory and kinematics. These kinematics are then used to investigate the 
effects of drag on the CME. We present the reconstructed CME (apex) kinematics for three 
events, one acceleration, one decelerating, and one with constant velocity. In Section 2 we 
describe the observations, data reduction, and the reconstruction and fitting technique. Sec- 
tion 3 includes a discussion of each event in detail and presents the reconstructed kinematics 
themselves. The implications of our results and our final conclusions are given in Section 4. 

2. Observations and Data Analysis 
2.1. Observations 



The trajectories of three CMEs were reconstructed using observations from STEREO 
SECCHI. SECCHI consists of five telescopes, the Extreme Ultraviolet Imager (EUVI), the 
inner and outer coronagraphs (COR1 and COR2), and finally the Heliosphereic Imager 
(HI1 and HI2). COR1 images the corona from 1.4-4.0 R , while COR2 images the corona 
from 2.5-15 R . Both of the coronagraphs take sequences of three polarised images which 



can be combined to give total brightness (B) or polarised brightness (pB) images (Howard 



et al. 2008a; Thompson et al. 2003). The HI instrument is a combination of two refractive 



optical telescopes with multi-vein, multi-stage light rejection system which images the inner 
Heliosphere from 4-89 degrees (Eyles et al. 2008). HI1 images the inner Heliosphere from 
3.98-23.98° (degrees elongation) in white light with a cadence of 40 minutes while HI2 images 
the Heliosphere from 18.68-88.68° in white light with a cadence of 2 hours. 

The three CMEs considered here were observed during: 2007 October 8-13 (CME 1), 
2008 March 25-27 (CME 2), and 2008 April 9-12 (CME 3). The observations were reduced 
using secchLprep from the SolarSoft library (Freeland & Handy 1998). This consisted of 
debasing and flat-fielding for all images. The COR1 and COR2 images were also corrected 
for vignetting, exposure time, and an optical distortion. The COR1 observations had a 
model background subtracted to remove static coronal features. The HI instrument has 
no shutter, and as such, these observations needed additional corrections for smearing and 
pixel bleeding. The pointing of the HI observations were updated using known star positions 
within the filed-of-view dBrown et al.|2008| ). Standard running difference images were created 
from the COR1/2 observations while a specialised running difference technique was used to 
suppress the stars for the HI observations (Maloney et al. 2009). The relative drift, due to 
satellite motion, of the star field between two successive HI images is calculated and then the 
earlier image is shifted to account for this motion removing a large part of the background 
signal. Figure 1 shows reduced observations from the 2008 March 28 event where the CME is 
simultaneously observed in both COR1 and COR2 in both the Ahead and Behind spacecraft 
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but only in from the Ahead spacecraft in HI. 




Fig. 1. Sample images of the 2008 March 25 event showing the CME in COR1, COR2, 
HI1, and HI2 fields-of-view. Observed from STEREO-A (top row) and from STEREO-B 
(bottom row). Note the absence of a clear CME signature in the STEREO-B HI1 and HI2 
images. 



2.2. 3D Reconstruction 



Each event was observed in either the inner coronagraph (COR1) or outer coronagraph 
(COR2) simultaneously by both STEREO-A and STEREO-B. From these images, the CME 



apex was localised via tie-pointing (see Maloney et al. 2009 Figure 1, Inhester (2006)). The 
trajectory was then reconstructed by tracking it through a series of images. In all the events 
presented, the CME was only observed in HI by one spacecraft, so an additional constraint 
was required to localise the CME apex. We therefore assumed that the CME continued along 
the same path with respect to solar longitude, as it did in the the COR1/2 field-of-view (i.e. 



travelled radially; Maloney et al. 2009). Figure 2 shows the derived trajectory for the 2008 



March 25 event. Once the 3D trajectories were derived, we calculated the height and then 
took numerical derivatives with respect to time to obtain the velocity and acceleration. 
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Fig. 2. — Reconstructed apex trajectory of the 2008 March 25 CME. The CME is tracked 
through COR1, COR2, HI1 and HI2 fields-of-view. For the COR1 and COR2 observations 
the CMEs position in reconstructed using triangulation. In the HI field of view we assume 
the CME will continue along the same path as it did in COR1 and COR2 (traces in x-y 
plane) in order to determine the CMEs position in 3D space. The x-axis points towards 
the Earth, the y-axis is is perpendicular to x-axis and in the ecliptic plane, the z-axis is 
perpendicular to both. 

2.3. Kinematic Modeling 

The kinematics were fitted, via a least-squares method, with a parametric model for the 
drag (Eq. [2]). In order to test which form of drag is most suitable (linear or quadratic), we 
fitted (Eq. [2]) with 5 set to 1 and then separately with 5 equal to 2. The kinematics were 
only fit during the time interval we believe that drag is at play and the observations are 
accurate. There was evidence for an early acceleration phase not attributed to drag which 
was not fitted. Also, events which were tracked far into HI2 field-of-view where identification 
of the CME apex become ambiguous were excluded from fitting. 

A number of the model parameters can be fixed from the observations, such as the 
CME height and velocity. We assume that the CME tends to the SW speed, which was 
taken to be where which the velocity plateaus. The model parameters obtained from the 
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fitting were then compared with previous results from Vrsnak & Gopalswamy (2002). From 
this comparison, we infer which model best reproduces the kinematics and hence is the most 
appropriate. Both the fast and slow CMEs (CME 1 and CME 2) were analysed using this 
method. The intermediate CME (CME 3) was fitted with a constant acceleration model to 
show that there was no significant acceleration involved. 



3. Results 



CME 1 and CME 2 were fit in three ways with a, /3, and S all allowed to vary (black 
line), with fixed 5 of two (magenta line) and one (orange line). In both cases the free fitting 
returned values that were not comparable to previous studies ( |Vrsnak 2001). The different 
fit parameters for the events are given in Table 1. CME 3 was fit with a constant acceleration 
model (black line). 



3.1. CME 1 (2007 October 8-13) 

Figure 3(a)-(c) shows the kinematics for the accelerating CME. This CME was first 
observed at 15:05 UT on 2007 October 8 off the west limb and was found to be propagating at 
an angle of 56° from the Sun-Earth line. Figure 3(a) shows the height of the CME. Figure 3(b) 
shows the velocity profile which clearly shows the CME is undergoing acceleration, initial 
velocity of ~150 km s~ x and final velocity of ~450 km s" 1 . There may be two acceleration 
regimes, an early increased acceleration phase (before 18:00 UT on the October 8) followed 
by a drag acceleration. The early acceleration can be attributed to a magnetic driving force 
and so was not fitted with the drag model. Later, when the CME reached the centre of 
the HI2 field-of-view, determining the front position becomes difficult so this region was not 
fitted. Figure 3(c) shows the acceleration profile of the event. The 5 = 2 (orange) fit givies 
the lowest chi-squared value. 



3.2. CME 2 (2008 March 25-27) 

The kinematics from the decelerating CME are shown in Figure 3(d)-(f). This CME was 
first observed at 18:55 UT on 2008 March 25 off the east limb and was found to be propagating 
at an angle of -82° from the Sun-Earth line. Figure 3(d) and (e) show the height and velocity 
profiles, the velocity profile clearly demonstrates the CME is undergoing deceleration. The 
CME had an initial velocity (HI1) ~800kms _1 and final velocity of ~375kms _1 . Due to 
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Fig. 3. — Left panel: Kinematics for CME 1, (a) height, (b) velocity, and (c) acceleration. 
Right panel: Kinematics for CME 2, (d) height, (e) velocity, and (f) acceleration. Vertical 
dashed line indicates separation between early and late phase acceleration. Horizontal dot- 
dash line indicates the inferred SW velocity. 
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the high speed of this CME, it was only observed in a small number of frames in C0R1 and 
C0R2. As a result, the kinematics were difficult to quantify in these instruments. However, 
there appears to have been an early acceleration feature. The deceleration in the HI1 and 
HI2 field-of-view continued until the CME reaches a near-constant velocity, and travels at 
this velocity throughout the rest of the field-of-view. Figure 3(d) shows the acceleration 
profile of the event. The 5 = 1 (magenta) fit gives the lowest chi-squared value. 



3.3. CME 3 (2008 April 9-12) 

In Figure 4 we show the kinematics of the constant velocity CME. This CME was first 
observed at 15:05 UT on 2008 April 09 off the east limb and was found to be propagating at an 
angle of -73° from the Sun-Earth line. Figure 4(a) shows the height of the CME Figure 4(b) 
shows the velocity profile which has a scatter about ~ 300 km s" 1 . Again, there may be some 
evidence in the COR1/2 observations for an early acceleration phase but due the events 
poorly observable features, at this early stage, it is hard to quantify this. The departure 
from the fit after April 12 20:00 UT is thought to be due to error in the reconstruction as the 
CME apex becomes to faint to identify. As this event shows no obvious acceleration it was not 
fitted with the drag model but with a constant acceleration model h(t) = h + v$t + l/2at 2 
(thin black line). Figure 4(c) shows the acceleration profile, the fit values (ho = 22 Rs un , 
t>o = 334kms" 1 , and a=— 0.18ms" 2 ) are consistent with no acceleration throughout the 
field-of-view. 

Table 1. Fit parameters for the accelerating events. CME 1, v sw = 450 km s -1 , 
v cme = 233 kms" 1 , h = 12 R . CME 2, v sw = 325 kms -1 , v cme = 702kms~ 1 , h = 44R . 



/3 s 



CME 1 (2007 Oct 8) 

Linear (magenta) 1.61e-5 -0.5 1.0 8.27 

Quadratic (orange) 1.28c-7 -0.5 2.0 6.74 

CME 2 (2008 Mar 25) 

Linear (magenta) 1.02e-4 -0.5 1.0 3.71 

Quadratic (orange) 6.38e-7 -0.5 2.0 17.63 
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Fig. 4. — Kinematics for CME 3 (a) height (b) velocity (c) acceleration. This event shows 
an early acceleration (to left if dashed line) which levels of to a scatter about typical solar 
wind speeds. This event was fit with a constant acceleration the resulting fit parameters 
are h = 22Rsun, t>o = 334kms™ 1 , and a = — 0.18ms -2 . The assumed solar wind value is 
indicated by the horizontal dot-dash line 
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Discussion and Conclusions 



We have shown it is possible to derive the 3D kinematics of features, the CME apex in 
this case, in the inner Heliosphere (~2-25OR ) using STEREO observations. The 3D kine- 
matics are free from the projection effects of traditional 2D kinematics but may contain ar- 
tifact from the 3D reconstruction method (e.g, Maloney et al.||2009 ) and other sources. Both 
of the accelerating events showed two regimes in the velocity profile, a low down (< 15 R Q ) 
early rapid acceleration (in comparison to later values), followed by a gradual acceleration 
far from the Sun (> 30 R ). The early acceleration is thought to be due to a magnetic driv- 
ing force, as the solar wind velocity low in the corona (v sw (< 10 R ) <268kms _1 , Sheeley 



et al. 1997) is lower than the velocity already attained by the CMEs in both cases. Here we 



assume that the later acceleration is due the interaction between the SW and the CME, as 
in each case the CME attains a final velocity close to typical values for the solar wind. 

Considering CME 2 in Figure 3(d)-(f), it can clearly be seen that the velocity levels off 
to a constant value typical of the solar wind. We interpret this as the CME reaching the 
local solar wind speed, as a result the force acting on the CME goes to zero. For CME 1 
Figure 3(a)-(c), the velocity initially increases, however, there is a plateau towards the end 
after April 11 6:00 UT which occurs at SW like speeds. The height measurements towards 
the end is very scattered and shows rapid increase. This is most likely due to losing the 
front to the background noise and triangulating a different feature. CME 3 propagates at 
a roughly constant velocity, which is consistent with the drag interpretation. The CME 
appears to have already attained the local SW speed and therefore is not accelerated. The 
fitting results show that a linear dependence produces a better fit for the fast event (CME 
2), while a quadratic dependence better fits the slow event (CME 1). The differing range of 
the interaction CME 1 ~12OR and CME 2~ 80 R Q may be explained by the suggestion 



that wide low mass CMEs are more affected by drag than narrow massive CMEs (Vrsnak 



et al. 2010) 



Reiner et al. (2003) suggest that for fast events, a linear model of drag better reproduces 



the kinematics, which agrees with our findings. Vrsnak (2001) also suggested that a linear 



dependence might be appropriate, however the quadratic form has been studied much more. 
From a theoretical perspective, a quadratic dependence corresponds to aerodynamic drag, 
while a linear dependence suggests Stokes' or creeping drag. It is not currently clear which 
model is more physically correct. The fit parameters obtained do not agree with those found 



by Vrsnak (2001) and while our values are not unphysical, it is not clear why they differ so 



much from the previous studies. 

The mechanism behind the apparent differing forms of drag, linear (5 = 1) and quadratic 
(5 = 2), for the slow and fast event are unclear. The application of any hydrodynamic theory 
to a CME, such as drag, may be missing vital physics. Could the magnetic properties play 
a role modifying the form of the drag (reconnection, suppression of turbulence, wave energy 



transport)? For example Cargill et al.| (1996) showed that depending on the orientation of 
the flux rope and background magnetic field (aligned or non-aligned) the drag coefficient 
can vary between zero and 3. They also found that the magnetic field of the flux rope is 
important in order for its survival as it propagates. Further which form of drag is correct 
for a CME in the SW, the low Reynolds number viscous dominated Stokes' drag or the 
high Reynolds number turbulence dominated aerodynamic drag? In order to address these 
questions a larger sample study is needed in order to verify these effects are recurring and 
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observable phenomena and also to build up the statistics. 

We have shown it is possible to derive the true 3D kinematics for a number of CMEs 
in the inner Heliosphere. Based on this we have been able to conclusively show that CMEs 
undergo acceleration in the inner Heliosphere, more specifically, that due to its range and 
strength this acceleration is believed to be the result of some form of drag. This drag 
acceleration has important implications for space weather predictions and for the analysis 
techniques which assume CMEs travel at constant velocity through the Heliosphere. The 
HI observations of CMEs in the Heliosphere provide a unique and limited opportunity to 
study the propagation of CMEs and to understand the coupling between the solar wind and 
CMEs. 

This work is supported by Science Foundation Ireland Grant No. 07-RFP-PHYF399. 
We would also like to thanks the STEREO/SECCHI consortium for providing open access 
to their data and technical support. 
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